Image-potential band-gap narrowing at a metal/semiconductor interface 

Ryotaro Arita, Yoshiaki Tanida 1 , Kazuhiko Kuroki 2 , and Hideo Aoki 
Department of Physics, University of Tokyo, Hongo, Tokyo 113-0033, Japan 
1 Fujitsu Laboratories Ltd., 10-1, Morinosato-Wakamiya, Atsugi, Kanagawa 243-0197, Japan 
2 Department of Applied Physics and Chemistry, University of Electro-Communications, Chofu, Tokyo 182-8585, Japan 

(February 1, 2008) 

GW approximation is used to systematically revisit the image-potential band-gap narrowing at 
metal/semiconductor interfaces proposed by Inkson in the 1970's. Here we have questioned how the 
narrowing as calculated from quasi-particle energy spectra for the jellium/Si interface depends on 
r s of the jellium. The gap narrowing is found to only weakly depend on r s (i.e., narrowing ~ 0.3 
eV even for a large r s = 6). Hence we can turn to smaller polarizability in the semiconductor side 
as an important factor in looking for larger narrowing. 

PACS numbers: 73.30.+y, 73.40.Sx 



I. INTRODUCTION 

Disrupted translational symmetry at surfaces and in- 
terfaces provides a potentially rich playing ground for 
many-body effects. A seminal proposal was in fact 
made by Inkson back in the 1970's, who proposed 
that a metal-insulator transition can take place around 
the interface.ETEl Classically, his argument is as follows. 
When a metal is placed on top of a semiconductor, an 
electron in the semiconductor feels an image potential 
which is the interaction between the particle and its im- 
age charge in the metal. This leads to a downward bend- 
ing of the conduction band bottom toward the interface, 
l/(4ez), where e is the dielectric constant of the semicon- 
ductor and z is the distance from the interface. Similarly, 
the valence band top is bent upward by the same amount. 

Quantum mechanically, this becomes as followso. The 
contribution of the correlation term to the self-energy of 
the electron in the semiconductor is similar (~ — l/(Asz)) 
for the conduction and valence bands. On the other 
hand, while the screened exchange term almost vanishes 
for the conduction band, it amounts to (~ l/(2erz)) for 
the valence band. As a result, while the conduction 
band bends downward like — l/(4ez), the valence band 
bends upward like l/(4ez). This kind of band bending 
can occur over short distances (z ~ O(10) A), while the 
usual Schottky barrier occurs over much larger distances 
(z ~ 0(100- 1000) A). 

After the proposal of Inkson was made, vari- 
ous studies for the band gap reduction or closure 
at metal/semiconductor (insulator) interfaces have been 
performed theoretically and experimentlvETEl Recently, 
Murata et aa have studied Ru(0001) /AI2O3 and have o 
served the band gap narrowing of AI2O3. Kiguchi et a, 
have studied LiCl films on Cu(001), and have found that 
3p level of CI shift up to the Fermi level as the number 
of LiCl layers decreases. 

As for the first ^principles many body calculation, 
Charlesworth et aln have calculated the quasiparticle 
electronic structure of Al/GaAs(110), and have shown 



the band gap narrowing for the first time. The amount 
of the band gap reduction turned out to be about 0.4 eV. 

A big issue remains however: which combination of 
metal and semiconductor will favor the local metal- 
insulator transition? As a first step toward such studies, 
we should investigate how the band gap reduction de- 
pends on the density of electrons (as represented by r s ) 
in the metallic side, which governs the electron correla- 
tion in that side. The value of r s in metals in fact extends 
over a wide range, 1.8 ~ 5.6, and the image-potential ef- 
fect in metals with greater r s is expected to be smaller 
than that for a smaller r s = 2.1, a value corresponding 
to Al and assumed by Charlesworthp. If the band gap 
narrowing still occurs significantly for metals with larger 
r s , we can turn our attention to the semiconductor side 
in optimizing the local metal-insulator transition. 

This is exactly the purpose of the present paper, i.e., 
to discuss the r s dependence of the band gap reduction 
quantitively. Since r s governs the dielectric response, the 
image-potential effect may well depend sensitively on r s , 
which is why we have to look into the dependence from 
first principles. For that purpose we have to go beyond 
the usual local density functional approximation (LDA), 
since we are talking about the effects of screening. So 
here we adopt what is called the GW approximation, 
which is roughly the RPA (random-phase approximation) 
+ LDA. 



II. MODEL 

To focus on the problem described above, we can sim- 
plify the metallic side into the jellium model. On the 
other hand, we have to have an atomistic model for the 
semiconductor side, since we are questioning effects oc- 
curring on the length scale of few atomic spacings on this 
side. So the model is depicted in Fig.l. To facilitate the 
band calculation, we adopt a repeated-slab model (pe- 
riodic boundary condition), in which the semiconductor 
slabs alternate with the jellium slabs. We calculate the 



band structure in the k x -k y space. 

When a semiconductor and a metal are put together, 
the Fermi energies in the metal and semiconductor have 
to be aligned in thermal equilibrium, which implies that 
some charges should flow across the interface resulting in 
a charged region on the semiconductor side in general (in 
the absence of impurity levels in the semiconductor gap) . 
This charge redistribution creates an electrostatic poten- 
tial, and this bends the valence and conduction bands 
of the semiconductor, which becomes significant at dis- 
tances z ~ 0(100 ~ 1000) A. Since what we want to look 
at is the physics on the scale of O(10 ~ 100) A, we can 
neglect this bending (unless the chargei-rearrangement is 
drastic). In fact, Charlesworth et alE have estimated 
the effect of the charge redistribution for Al/GaAs(110), 
where electrons flow from Al with a higher Fermi energy 
to the lower conduction band bottom in GaAs, and found 
the effect to be negligible. However, when the charge re- 
arrangement occurs drastically, this one-body effect can 
become serious. Thus, to focus on the image-potential 
effect we have to exclude the one-body potential effect 
carefully. For this purpose, we set the Fermi energy of the 
metal inside the energy gap of the semiconductor(Fig.2). 
As we shall show, such a situation is indeed realized if we 
take an appropriate value of r s for the jcllium model. 

As for the semiconductor side, we have employed a 
slab, which consists of five layers of Si atoms stacked in 
the [001] direction. To get rid of complications arising 
from dangling bonds, the edges are terminated by hydro- 
gen atoms. For the structure of the hydrogen-terminated 
surface, we assume for simplicity a non-reconstructed 
one (p(l x 1)). With this assumption we have opti- 
mized the structure imposing the mirror-plane symme- 
tries (along [100] and [010]) and an inversion+mirror 
symmetry ([001]). The hydrogen atoms are allowed to 
relax in any directions. The size of a supercell is 7.1 8 2 x 
30.48 a.u. 3 along the (x,y, z), with the thickness of the 
jcllium being 12.51 a.u. 

The band gap narrowing or closure in the semiconduc- 
tor is probed here by identifying the character of the wave 
functions for various bands for the repeated slab model: 
by concentrating on the bands whose wave functions have 
their amplitudes primarily on the semiconducting side we 
can define the gap of the semiconductor. 



III. METHOD 

Band-structure calculations are usually performed 
within the framework of the LDA. In this formal- 
ism, many-body effects are represented by the so-called 
exchange-correlation potential, which is a functional of 
the electron density. In practice, this potential is approx- 
imated as a function of the local density, and we have the 
LDA. 

While the density-functiow&l formalism is shown to be 
rigorous for the ground statet3 and LDA gives reliable in- 



formation about the ground state properties for various 
electron systems, it is well-known that these approaches 
are not useful for excited states. In fact, LDA usually 
underestimates the band gap of semiconductors and in- 
sulators. Moreover, LDA cannot be applied to the cases 
where the electron density varies in space. 

Still, LDA wave functions are usually good approxima- 
tion to quasi-particle wave functions!!! Since excitations 
can be described by many-body perturbation theory, it 
should be quite a good starting point to adopt the LDA 
wave functions as the basis for the many-body perturba- 
tion theory in determining the self-energy and spectrum 
of the quasi-particles. 

For the calculation of the self-energy, various approxi- 
mations have beeHrjdeveloped. Among them, Hedin's GW 
approximation! 12 !] 13 ! often gives excellent quasi-particle 
energies in bulk . sem iconductors with a comparatively 
simple formalism! 1 HI 14 !. The GW approximation essen- 
tially amounts to the RPA in the LDA formalism, so we 
have adopted this method to study the image-potential 
band-gap reduction. 

A. LDA 

So the first task in the present study is to perform 
an LDA calculation to obtain the eigenwavefunctions for 
the system described above. We adopt the exchange- 
correlation functional introduced by Perdew and WangE_3, 
and eigenwavefunctions are expanded by plane waves up 
to a cut-off energy of 16Ry. As for the atomic pseudo- 
potentials, soft«-*orm-conserving pseudo-potentials in a 
separable forrrJl3 are employed. The atomic configu- 
rations and the corresponding electronic states in the 
groundj-states are obtained with the conjugate gradient 
schemeEZI. 



B. Self-energy correction in the GWA 

We then proceed to the GW approximation 
(GWA)calculation. The central idea of GWA is to ap- 
proximate the self-energy operator £ by 

£(r,r';u;) = — / du'G{r, r'; u + w')W(r, r'; uj')e luJ ' s , 
2ir J 

where S is an infinitesimal positive time and W is the 
screened Coulomb interaction, 

W{v,v'-lo) = [ dv" - 1 

J e(r , r , w)\r — r" 

where e is the dielectrkiiunction. 

Recently, Rojas et a$3 proposed a new implementation 
of GWA, the space-time approach, which is described in 
detail by Rieger et a££l. The quasi-particle calculations 
can be performed either in reciprocal space as a function 
of the frequency, or in real space as a function of the 
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imaginary time, and the central idea in this method is 
to choose the representation that minimizes the compu- 
tations required to evaluate the basic GWA quantities. 
This approach enables us to study larger systems. 

The actual computational steps in this method are as 
follows. First, split the self-energy into a bare exchange 
part Yi X and an energy-dependent correlation contribu- 
tion T, C (E). The former can be calculated from 
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Here G is the reciprocal vector, and 



M^(k,q) = J $ t , k _ q (r)e- i (i+ G ) r $ m:k (r)dr, 

where ^^^-q is the wave function in the valence band 
of the semiconductor, while $ m , k is the wave function in 
the m-th band. 

To evaluate the energy-dependent self-energy T, C (E), 
we first construct Green's function in real space and 
imaginary time, 
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where and e n are LDA wavefunctions and eigenvalues. 
Due to a rapid decay of the exponentials, the convergence 
against the cut-off in X)nk° CC * s mucn better than that in 
the real frequency formalism. 

Next, the RPA irreducible polarizability x°( r i r '; * T ) ls 
calculated in real space and imaginary time, and Fourier- 
transformed to reciprocal space and imaginary energy- 
Then, the symmetrized Hermitian dielectric matrixc2l 
£g,g' (k, iuS) is constructed, and inverted for each k point. 

Then the screened Coulomb interaction is calculated 



W G g> (k, iuS) 
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and is Fourier-transformed to real space and imaginary 
time. The self-energy operator can be calculated as 

S c (r, r'; it) = iG LDA (r, r'; ir)W(r, r'; it). 

Finally, we evaluate the correlation contribution, 
£ (it) — (kn|S c |kn). After this is Fourier-transformed 
to imaginary energy, we perform analytic continuation 
onto the real energy axis with the Pade approximationEi 
In the present study, we consider a 6 x 6 x 24 grid for 
the unit cell, and a 6 x 6 grid for (k x , k y ). The time grid 
is spaced by 5t = 0.3 a.u. over the range of 13 a.u. We 
have taken up to 253 states to construct Green's function. 



IV. RESULT 



Let us move on to the results. The quasi-particle spec- 
trum e® p is obtained as 



e QP _ £ LDA 



V X L C DA + S x + E c . 



In the following, we discuss the effect of the metallic layer 
for each term in e® p . 



A. The LDA calculation 

We first show the band structure and the (squared ab- 
solute value of) the wave functions obtained with LDA 
for r s — oo, i.e., no jellium (Fig. 3), r s — 6 (Fig. 4) and 
r s = 4 (Fig. 5), respectively. In the absence of the metal 
(r s = oo), we can see that the valence band top is at T, 
while the conduction band bottom lies around K (and 
T). Hereafter, we focus on the energy shifts of the va- 
lence band top around T and the conduction band bot- 
tom around K caused by the close contact with the jel- 
lium. 

When in vacuum the LDA band gap across K and V 
is 1.93 eV. When we introduce the jellium with r s = 6 
or 4, the Fermi level still lies across K and T, i.e., only a 
small amount of electrons flow from the jellium into the 
Si side, so the effect of the charge rearrangement is al- 
most completely absent. For smaller values of r s , on the 
other hand, we can show that electrons in the jellium do 
flow into the Si side, and the electronic band structure 
changes drastically. Thus we focus here on the case of 
r s = 6 and 4. 

For r s = 6(4), the band gap across T and K reduces to 
1.85 (1-79) eV. The characters of the wave functions at 
these points are mainly Si and do not change, as we can 
see in Fig. 3, 4, and 5. In addition we notice that there is 
a state which emerges around T crossing the Fermi level 
for r s = 4, 6 (the bands represented by bold lines in the 
figures). If we examine the character of wave functions 
on these branches in Figs. 4,5, they reside well within 
the metallic side, so we exclude them from our argument 
on the band-gap narrowing in the semiconductor side. 



B. Vxc and the direct exchange term E x 

Now we are in position to analyse the band gap reduc- 
tion term by term. Here we consider Vxc and the direct 
exchange term E x , the energy independent part of the 
GW correction to the LDA eigenenergy. The expectation 
value of the exchange-correlation potential is obtained by 
(V£ DA ) = / |$ k ,„(r)| 2 V^ DA (r)dr. We have found that 
V xc (r) LDA in the Si region has similar values for r s = oo, 
6 and 4. Since the characters of the wave functions do 
not change as we have seen in Figs. 3, 4, and 5, we can 
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expect that (V^f A ) does not change significantly when 
we introduce the jellium between the Si system. 

On the other hand, the value of the matrix elements 
between the metal wave functions and the semiconductor 
wave functions, M'^Ck, q), which governs £ x , are small, 
since the states of jellium character and Si character are 
well separated in real space (see Figs. 4, 5). Therefore, we 
can expect that T, x does not change significantly when 
we introduce the jellium. 

In fact, -Vxc + S x at K is 3.99 -> 4.00 -► 3.98 in 
eV for r s — co — > 6 — ► 4, so the shift is negligible. At 
r, -V X c + S x is -2.09 -> -2.00 -» -1.93 in eV for 
r s = oo — ► 6 — > 4. 



C. Energy-dependent E 

Finally, let us discuss the correlation contribution, S c . 
In Fig. 6, we show the imaginary frequency dependence 
of E c at T and K. 

We have then performed an analytic contmuatiorupnto 
the real energy axis with the Pade approximation!^. In 
Fig. 7, we show the real frequency dependence of £ . We 
can see that there is an almost linear dependence on oj. 

When E c is linear for small u>, the GW correction to 
the LDA spectrum reduces to 

A = -^-{*r*\E c (eg A ) + - V XC \^) 

where 

Znk = i - -^($ nk |£ c H|$„ k )L =e ^A. 

From Fig. 7, we can see that Z nic = -0.2 ~ -0.3. At T, 
A is estimated to be -0.75 -> -0.60 -> -0.44 in eV for 
v s = oo — ► 6 — > 4. At K, A is 1.28 -> 1.15 -> 1.09 in eV. 
Thus the band-gap reduction due to the presence of the 
jellium amounts to as large as pa 0.3 eV for r s = 6 and 
pa 0.5 eV for r s = 4. 



V. DISCUSSIONS 

To summarize, we have studied, with the GW 
approximation and a character-resolved band anal- 
ysis, the image-potential band-gap narrowing at a 
metal/semiconductor interface by calculating quasi- 
particle energy spectrum of the jellium/Si interface. For 
the values of r s = 4 — 6 of jellium studied here the elec- 
trons or holes do not flow from the metallic to semicon- 
ducting side, i.e., the Fermi energy of the jellium lies 
within the energy gap of Si, so the one-body effect due 
to the charge redistribution is absent. We have found 
that a significant band gap narrowing of pa 0.3 eV occurs 
for r s as large as 6. 



So we can concentrate on the semiconducting side 
to realize larger gap-narrowing effects or a local metal- 
insulator transition. If the dielectric constant of the semi- 
conductor is small, the image potential effect will become 
stronger, so that we may expect larger band gap narrow- 
ing. However, the system with a small dielectric function 
usually occurs in materials with large band gap, so that 
the realization of the band-gap closure becomes a trade- 
off. Furthermore, the energy gap of the semiconducting 
layer may depend on the surface structure (e.g., whether 
it is terminated by hydrogen atoms or dangling bonds 
form dimers, etc.), and their effect is also non-trivial. 
These are interesting future problems. 

Finally, while the band gap narrowing or closure is of 
a fundamental interest in its own right, we can further 
raise a very strong motivation for considering a metaL 
lized semiconductor surface. There is a long historyEj 
of proposals for superconductivity in conducting systems 
in close contact with polarizable media. Littleo prp-i 
posed this for one-dimensional systems, then Ginzburgcil 
extended this to two-dimensional systems. Allender, 
Brey and Bardeenn3 studied this in detail, with metal- 
semiconductor structures in mind, which .-was subse- 
quently criticized by Inkson and AndersonEj. A sum- 
mary of the situation by Zharkovn3 is that, while the 
criticism is correct for the usual form for the dielectric 
function, the situation may be resurrected for unusual 
dielectric functions. The background to all this is that 
for the polarizable-medium mediated superconductivity, 
the conduction layer should be very strongly coupled to 
the semiconducting layer, ideally with strong chemical 
bonds such as covalent ones. So the metallized semicon- 
ductor surface, with the band-closure mechanism, should 
be one ideal realization of this. 
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FIG. 1. The model studied in the present work. The slabs, 
each of which consists of Si layers terminated by H atoms, 
sandwich the jellium (metallic) region. 
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FIG. 2. Two cases where there is a charge transfer across 
the metal-semiconductor interface (top panel) or there is no 
charge redistributions with the Fermi energy of the metal ly- 
ing within the gap of the semiconductor (bottom). 





FIG. 3. The band structure of five Si layers terminated by 
hydrogen atoms in vacuum. The squared absolute value of 
the LDA wavefunctions at the valence top and conduction 
bottom are also shown. 




FIG. 4. A plot similar to Fig. 3, when the jellium with 
r s = 6 is attached to the Si layer. The branch whose am- 
plitude is localized in the jellium region is denoted by dotted 
lines. 
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